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Abstract 

A brief overview is given of some recent advances in charged-composite particle 
scattering. On the theoretical side, I address the three-charged particle wave func- 
tion asymptotics, the nonperturbative investigation of the long-range behaviour of 
the optical potential, and the question of the compactness of the kernels of the 
momentum space integral equations for three charged particles. Among the more 
practical developments, I report on results of numerical calculations of so-called 
"triangle" amplitudes, a new, simple and very efficient higher-energy approxima- 
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effects in proton-deuteron elastic scattering with realistic nuclear potentials. 



1 Theoretical Developments 



1.1 Asymptotic wave function for three free charged particles 



Knowledge of the asymptotic boundary condition on the three-free charged particle wave 



function is required not only when attempting to solve the Schrodinger equation above 
the ionization threshold, but also when investigating asymptotic properties of various 
interesting quantites like the optical potential or the behaviour of the kernel of momentum 
space integral equations. I, therefore, start by briefly recapitulating some important 
aspects of the asymptotic behaviour of the solutions of the Schrodinger equation. 

Consider three distinguishable particles with masses m v and charges e vi v = 1,2,3. 
I use Jacobi coordinates: k a (r a ) are relative momentum (coordinate) between particles 
(3 and 7, q Q (p a ) the relative momentum (coordinate) between particle a and the center 
of mass of the pair (fij). The on-shell momenta are denoted by (k Q ,q a ) so that the 
on-shell relation reads as E = q^/2M a + k^,/2fi a , with fi a = m / gm 7 /(m^ + m 7 ) and 
M a = m^mp+m^) /(rria+mp+rriy) being the appropriate reduced masses. Furthermore, 
V a = V£ + V£ is the short-range plus Coulomb potential acting between particles (3 and 
7, and fj a = r] a (k a ) = epe^pa/ka the corresponding Coulomb parameter. Additional 
notation: E + = E + iO, e a p = —tp a = +1 if (a,/3) = cyclic ordering of (1,2,3); finally, 
unit vectors are characterised by a hat: k = h/k. 

The asymptotic solutions of the Schrodinger equation for three asymptotically free 
particles in the various regions of configuration space are known: 

In fl : n, Tii r 3 — * 00 j but not r u /p u — > for v — 1, 2, 3 (Redmond, as cited in Ref. M): 



(r O ) ~ e i s -ot- r <x+<lc t -p oi ) TT e i7?„ln(ft„rv-k„-r„) r-y 

u=l 

In n 

a'- Pa — y 00, r a / Pa — y 0, for a — 1, 2 or 3 (0, with some refinements given in |j]) 
' ( - +) (r o ) « 
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Here, ip^ jp )( r a) is continuum solution of the two-body-like Schrodinger equation 

describing the relative motion of particles (3 and 7 with local energy E a (p a ) = k^ l (p a )/2fj JOl , 
where (e.g. A/3 = fj, a /m 7 , with (3, 7 7^ a) 

k a (p a ) = k a + , a a (p a ) = - 2^ VuK — . (4) 

P° 1 - e au p a • K 

Its parametric dependence on p a is a manifestation of long-ranged three-body correlations; 
thus it is, in fact, a three-body wave function, the influence of the third particle a however 
being confined to a shift of the relative momentum of particles (3 and 7 from its asymptotic 
(for p a — > 00, i.e. particle a is infinitely far apart) value k Q to the local value k a (p a ). 



1.2 Long-Range Behaviour of the Optical Potential in a Three- 
Body System 

Elastic scattering processes with (charged) composite particles of the type a + {(3, j) m — ■> 
a + (/3, j) m can formally be described by means of a single one-channel LS equation for 
the elastic scattering amplitude 

z - Q^/2M Q + E am 

where the plane wave matrix elements of the optical potential (OP) operator are given as 

Km.amCq'a; ^) = (q'J (VwlK + V a Q am G(z)Q am V a \lp am ) |q a ) 

static potential nonstatic part of OP 

The following notation is used: \ip am ) denotes the target bound state wave function to 
binding energy E am ; G(z) = (z - H - J2 V V v )~ l and G c (z) = (z - H - J2u K C ) _1 are 
the resolvents of the full and the pure Coulomb three-body Hamiltonian, and G a (z) = 
(z — H — Va)^ 1 the one of the channel Hamiltonian. P am = \ip am )(ip am \ projects onto 
the target state, and Q am = 1 — Pam onto the orthogonal complement. Furthermore, 
V a = Vcf + Vu = Ei/^« Vis is the channel interaction; Q Q is the momentum operator with 
eigenvalue q Q . 

The question of solvability of (|5|) depends on the singular behaviour of V°^ am in the 
limit that the momentum transfer A a = q' Q — q Q goes to zero. In leading order, the 
latter is caused by the Coulombic part V£ of the channel interaction, i.e., it does not 
depend on either the short-range part or on the internal interaction V a . Of course, a 
certain behaviour of the optical potential for A a — > implies a corresponding asymptotic 
behaviour in coordinate space. 
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As is well known, the static potential has, in the limit A a — ► 0, only the trivial 
Coulomb-type singular behaviour, which for a spherically symmetric target looks like: 

ystat ii \ Ao-+o 47reQ.(e/3 + e 7 ) -ystat / \ pq^oo e a (e/3 + e 7 ) , . 

'am,am\4ai "a/ a 2 'ara,nm \rixl ■ V / 

The behaviour of the nonstatic OP has been known for a long time, but only in 2 nd 
order perturbation theory which in the present language is equivalent to approximating 
in ((^) the three-body (G) by the channel resolvent (G a ), or in adiabatic approximation. 
On the energy shell and below the ionisation threshold one has (with A a = — q Q ) 



^«m,am\Ha)4ai 1J +^ > r y am,an 



ara,am\raJ n 4 ' V / 

Here, a is the so-called static dipole polarisability of the composite particle. Two questions 
arise immediately: 

(i) Does the fundamental result (||) hold also for the exact nonstatic OP, i.e., even after 
all terms of the perturbation expansion of G(z) are summed up? (For E < 0, this has 
been answered in the affirmative, though not fully rigorously, in ||.) 

(ii) Does it hold also for E > 0? (Within perturbative approaches the answer was yes, 
but with an energy- dependent "a" in ||, and with the standard a in []7|.) 

The behaviour of am for A a — > has been investigated in || nonperturbatively 
and for all energies, by inserting the spectral decomposition of G C (E + ) with both two- 
and three-body intermediate states; for the latter the asymptotic three-charged particle 
wave function in Q a has been used. Off the energy shell we find: 

V£U(<& q«; E ) A =° Ci A * + o (A Q ), (9) 

with 
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D k o m | 2 + |A Q -D 
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Cl 4 [£ + - q l/2M a - E an ] + J (2vr)3 [ E+ - q^/2M a - kf/2^ ' ' (10) 

and (N a = e a/3 e a /i a [e 7 /m 7 - ep/mp]) 

D nm = N a J dr a ijj* an (r a )r a ip am (r a ), D k o m = N a J dr a ^J*(r a )r a ip am (r a ). (11) 

On the energy shell, am is seen to depend on the momenta only via A a , and no longer 
on the energy. Hence, in coordinate space it is a local, energy-independent potential 

V^JpJ^-A + of^V (12) 

Z Pa \PaJ 

with the polarisability as known from the perturbative approaches: 

q = 2E r ^• P " m|2 t Wt^ T J^"^/ t (13) 



dk° 
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+ kf/2ij a 



That is, no "renormalisation" of a arises from the higher order terms in the perturbation 
expansion of G(z), and all dependence on E which is present in the off-shell "strength 
factor" Ci(E) has disappeared. Note that along the same lines the existence of a new 
nonrelativistic contribution to V°^ am ~ 1/p 5 has been established recently in ||. 
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1.3 Compactness properties of the kernels of momentum space 
integral equations for three charged particles 

Three-body integral equations of the Faddeev type can not be used for Coulomb-like 
potentials, because of the occurrence of singularities in their kernels which destroy the 
compactness properties known to exist for short-range interactions. Up to now, only for 
energies below the breakup threshold had attempts been successful to obtain integral 
equations with compact kernels, by singling out from the original kernel the so-called 
two-particle Coulomb singularity, in a form such that it could be inverted explicitly |T0| . 

To investigate the behaviour of the kernels for positive energies, we use the rigorously 
equivalent formulation in terms of an effective-two-body theory flll| . Here, the transition 



amplitudes for all binary processes satisfy a closed set of coupled LS-type equations 

3 

Tgn,am{z) = Vffn.amjz) + / J Vffn.w (z) Go\v,rs ( z ) %/s,am (z) ■ (14) 

v=\ r,s 

Without loss of generality, the short-range potentials can be taken as separable potentials: 
^cf = Em \Xam)^am(Xam\- The effective potential matrix elements are given as 

V/3n,am(q!f3, <la] z) = (cfp, Xf3n\ G° (z) ~ 5p a G^(z) \Xam, q«). (15) 

The effective propagator matrix has simple poles in its diagonal for those m which corre- 
spond to bound states: ^o;«,mm(qai z) ~ l/(z — q 2 a j2M a — E am ). Note that as a result of 
our choice of the form of V£ , these expressions contain only pure Coulombic quantities. 

Define D (z) = z - k^/2fjtp - qf/2M p = z - k 2 a /2fi a - q 2 a /2M a , with k' p = \q a + 
[ipqUrrijl, k a = |q^ + // a q Q /m 7 |, r]p = T]p{k'p) and r\ a = T] a (k a ). Then, we have now 
proved |T2j that even for positive energies, provided all three particles have charges of 
equal sign, the leading singularity of the nondiagonal off-shell effective potential is a 
branch point at D (z) = 0: 

Vto^, q«; z) Do( ^° l/D (z)^ + ^\ for (3 ^ a. (16) 



This compares with the simple pole occurring in the 'pole' amplitude V|°^ m (q^, q, 
<Wq/3,X/?J Go{z) \Xam,<la) ~ l/D (z). On the energy shell, it goes over into 

Vpn, am {^ q.; E+) 1/D^^*\ for (3 ± a, (17) 
where we have defined the bound state Coulomb parameter via = epe~ffi a / n am , with 



\z\ 



1 - 2fi a \E am \, and similarly for , and denoted by D = k 2 t /2fi a + \ E am \ the on-shell 
restriction of Dq(z). Since neither the pole nor the branch point singularity can coincide, 
for real values of the momenta, with the effective propagator pole, both are harmless. 

Preliminary results indicate that the diagonal kernels develop on the energy shell 
a nonintegrable singularity which is, however, of same two-body type (" center-of-mass 
Coulomb sing ularity") as that found in [JTOj for E < 0, and in [0, H (see also ||15|| ) for all 



am 



energies within the screening approach; thus, it can be treated by explicit inversion. Apart 
from that the next strongest singularity occurs off shell and turns out to be 0(A~ 5//2 ) for 
A Q — > 0. Thus, it is integrable. 

Hence, it appears that after a few iterations the kernels are compact even above the 
breakup threshold provided the charges of all particles are of equal sign. In contrast, if 
some of the charges have opposite signs the kernels develop nonintegrable singularities 
which can destroy the compactness properties. 
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2 Practical Developments 



2.1 Energetic Collisions of Charged Projectiles with Atomic 
Bound States 

Scattering of an elementary charged projectile off a two-charged particle bound state such 
as hydrogen atoms, positronium, etc., at higher energies can be described by the first few 
terms of the multiple-scattering expansion of the three-body transition operator: 

^,am(^. Qo) = (q|s> fpptJl {ha G o(E+) + Sp v 6 va T^(E + ) H } \ip am , q a ). (18) 

Here, G (E + ) = (E + iO — Hq)' 1 , the two-body Coulomb T-operator acting between 
particles a and (3, and 5p a = 1 — 5p a . The momenta satisfy the on-shell condition E = 
q 2 j2M a + E am = q'$/2M p + E 0n . 

The first ('pole') term of ([Tj|) describes the elementary one-particle transfer. The 
first-order rescattering ('triangle') amplitudes ($ = Z(q«,q a )) 

£ M«<* I (^n I If (£ + iO) | ^ Qm > | q Q > = { ^^"3^0 1 P -t °L _ (19) 

contribute to both direct ((3 = a) and exchange (/3 7^ a) scattering. They (as well as all 
higher-order contributions) contain T c describing the rescattering of projectile a off each 
charged target particle. Hence, their calculation is difficult and time-consuming (this 
is particularly so if the rescattering particles have charges of opposite sign, because of 
the occurrence of the infinity of bound states). Thus, usually the approximation — > 
is made, yielding the so-called Coulomb-Born Approximation (CBA), to be denoted 
by M-^ nam and M.„ nm - But for atomic reactions the CBA is known to fail badly (for 
attractive and repulsive cases), except for E — > 00 (see |16[ and references therein). 

We have succeeded to numerically calculate all (direct and exchange) triangle ampli- 
tudes for all E and $ for arbitrary wave functions, for repulsive and attractive intermediate- 
state rescattering (for the latter, a 'new' representation of the attractive Coulomb T- 
matrix for E < was developed which has the bound state poles displayed explicitly in a 
simple, numerically convenient manner). Many concrete results are discussed in [17]]. 

Moreover, we have derived a new approximation in the form of a 'renormalized' CBA: 

M^Jfi, E) « E) A<L(tf, E), 7 ^ a, (20) 

A< am (tf, E) « n ( ;l am ($, E) M^ >am (#, E), /3? a , (21) 

<T?(s) (fl TP) _ i^ m \ ) / / 22 n 

'S^m^'^J „, (d) ( d)' 7f «, l^J 

A ^ 7 A 2 A 2 -\- ( k 4- «■ ")2 ~ 2M h 

2) = ^ 2 . w , 77^^«. (23) 

A2 A 2 , , „ ^2r 2^ ^ 



with 



Here, the following notation is used: = q^ — X a q a /Xp, X a = p, a /mp, Xp = fip/m a ; 

n 2 am = 2fi a \E am \ and similarly for K 2 pn . Finally, r]^ e) = e a ep^ / {2^ 1 (E + -k 2 {de) /2M 1 )} 1 / 2 . 
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The quantities k^(E), k^(E) are real, A(E), B(E) are real for rj^^ real, and 11(E) is 
complex. All of them are independent of and are explicitly given in terms of simple 
functions. We note in parentheses that we have also derived approximate analytical 
expressions for Ai v (for arbitrary bound state wave functions). When inserted into (]2lD 
they yield approximations for the triangle amplitudes containing no quadratures at all 
(however, their range of validity is somewhat limited). Because of their simple structure, 
both types are very easy to use for theoretical as well as numerical purposes. 

Selected applications: Consider E > k 2 de J2M y so that r]!^'^ are real. Then (M.,71 
without channel indices refer to both 'diagonal' and 'nondiagonal' quantities): 



1. The whole ^-dependence of 



M 



rpC 



is given by that of 



M T °($,E) 



1Z is) (E) 



M TC ($, E) 



(24) 



the CBA fails with respect to the magnitude only, insofar as \1Z^(E)\ differs from one. 
2. With the penetration factor Cq = 2nr]^ d ' e ^ /[exp{2TiT]^ d ' e ^} — 1] pertaining to the 
intermediate-state Coulomb scattering one finds 



M 



rpC 



(d,e) 



cg(i + ofo 



(d,eY 
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Vm, n, 



(25) 



Since Cq < 1(> 1) for e a ep > 0(< 0), ( |25|) not only quantifies the charge-sensitivity of 



M. T ° \ as compared to the insensitivity of lA^^"" | (i.e., Ai v ^ (e a ep > 0) 



-M 



e a e/3 < 



0)), but it also provides a simple method to quantitatively estimate the former. 
3. Since C 2 (e Q e^ > 0) = CQ 2 (e a e p < 0) (l + 0(r]^ 2 )) one derives 



M^ nm (e a e p > 0) 



M^ nm (e a e p < 0) 



(26) 



relating the 'direct' triangle amplitudes for processes with opposite signs of the charges 
of the rescattering particles. 

4. The dependence of fc(d, e ) o n n> and m becomes negligible for sufficiently large E, render- 
ing 7]^'^ and Cq practically independent of all bound state characteristics. Consequently, 
relation (|5|), and all results derived from it, become universal, i.e., state-independent. 

Many numerical tests of the approximate triangle amplitudes have been performed 
with e ± ,p,p as projectiles and H, Ps, (p,p) as targets Quite generally one finds for 



both direct and exchange reactions that they reproduce the exact triangle amplitudes 
already at 1 keV incident energy for light, and at 100 keV for the heavy projectiles, to 
within a few percent, for practically all scattering angles, the agreement becoming even 
better with increasing energy. This is to be contrasted with the CBA's: not only are they 
real, but their magnitudes for light projectiles at 1 keV are off by ~ 40 - 100 per cent, 
the improvement with increasing energy being slow only. Similarly for heavy projectiles. 



2.2 Proton-Deuteron Scattering with Realistic Potentials 

Ever since the development of the exact few-body theories, the investigation of the 
nucleon-deuteron (Nd) elastic scattering and breakup reaction has been at the center of 
interest, as their principal field of application. For processes with neutrons as projectiles 
(nd), highly sophisticated calculations are available nowadays, using realistic nuclear po- 
tentials including three-body forces. Apart from a few noticeable and as yet unexplained 
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failures, they provide a good to very good description of all the available observables (for 
a recent review see [fL8|| ). However, experimental nd data are rather sparse and are lacking 
the desired accuracy. Hence, usually nd calculations are compared with pd data which 
makes too good an agreement rather questionable. 

For proton-induced reactions (pd), on the other hand, data are abundant and of excel- 
lent precision, but similarly advanced calculations have hitherto been restricted to E < 
|H| pp| (calculations with simple nuclear potentials, which have been performed for all 
energies, give in general "only" semiquantitative agreement but they are usually very 
accurate in explaining observed differences between pd and nd data, see [F2T], and 
references therein). The reason for this lack of sophisticated pd calculations above the 
breakup threshold is obvious: In coordinate space, imposing Coulomb boundary condi- 
tions in the whole configuration space is a very difficult task. In momentum space, though 
the screening and renormalisation approach [13|, [14], [15| is in principle straightforward, its 
application is very computer time consuming (for a pedagogical introduction to the theory 
and a rather complete list of references, both on theory and calculations, see [|15|1). 



We have now succeeded to obtain for the first time results for cross section and po- 
larisation observables above the breakup threshold for a "realistic" (Paris) potential [[23 



within the screening and renormalisation approach. In fact, we use a separable represen- 
tation of the latter (so-called PEST1-6) known to provide an excellent approximation to 
the original local potential |24| . Up to now the nucleon-nucleon interaction has been taken 



into account in the states 3 Si — 3 Di, ^So and in all P waves. When comparison is made with 
selected experimental data at 5 and 10 MeV [^, p6| , we find good although not perfect 



agreement of the calculated with measured quantities. Part of the discrepancies (in the 
polarisation observables A y and iTu) exist already in the nd case, implying that there, as 
well as in our calculations, some aspects of the nuclear force are still missing. Other dis- 
crepancies may arise from the restriction in the number of nucleon-nucleon partial waves 
taken into account, or from a possible influence of the - hitherto neglected - three-body 
forces. However, the general trends, in particular concerning the relation between neutron 
and proton data, are well reproduced. Clearly, the observed agreement of nd-calculations 
with pd-data, in particular for iTu and also to a somewhat lesser extent for A y , is purely 
fortuitous and does not arise from smallness of Coulomb effects in that observable. 

A final remark concerns the fact that we have taken into account the Coulomb inter- 



action in CBA, which for the present case has been estimated in fll7| , jig] to be accurate 
to better than 1%. 
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